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The level surfaces of solutions to the eikonal equation define null or characteris- 
tic surfaces. In this note we study, in Minkowski space, properties of these surfaces. 
In particular we are interested both in the singularities of these "surfaces" (which 
can in general self-intersect and be only piece-wise smooth) and in the decomposi- 
tion of the null surfaces into a one parameter family of two-dimensional wavefronts 
which can also have self-intersections and singularities. 

We first review a beautiful method for constructing the general solution to the 
flat-space eikonal equation; it allows for solutions either from arbitrary Cauchy data 
or for time independent (stationary) solutions of the form S = t — Sq(x, y, z). We 
then apply this method to obtain global, asymptotically spherical, null surfaces 
that are associated with shearing ("bad") two-dimensional cuts of null infinity; the 
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surfaces are denned from the normal rays to the cut. This is followed by a study 
of the caustics and singularities of these surfaces and those of their associated 
wavefronts. We then treat the same set of issues from an alternative point of view, 
namely from Arnold's theory of generating families. This treatment allows one to 
deal (parametrically) with the regions of self-intersection and non-smoothness of 
the null surfaces, regions which are difficult to treat otherwise. 

Finally, we generalize the analysis of the singularities to the case of families of 
characterisitic surfaces. 
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I. Introduction 



The study of the propagation of electromagnetic wavefronts, in the high- 
frequency or geometric optics limit, is almost ubiquitous in physics; it is a 
basic staple of elementary physics courses, it arises in the practical area of 
optical equipment, in applied subjects, too numerable to mention in detail, 
involving materials with variable refractive index , in atmospheric and astro- 
physical studies. They have been a prime illustration of V.I. Arnold's theory 
of Lagrangian and Legendre maps.BBifl In a different guise, similar prob- 
lems arise in catastrophe theory. In addition to the various applications to 
more standard physics problems, they also play a most fundamental role in 
General Relativity, e.g., the continuous propagation of the two-dimensional 
wavefronts, i.e., the one parameter family of evolving wave fronts, form null 
(or characteristic) three-surfaces that are determined by the dynamics of the 
curved space-time in which the wavefronts propagate.!!! In this context it 
forms the basis for the theory of gravitational lensing.! The converse state- 
ment is also true, namely that sets of null surfaces define, up to a conformal 
factor, the space-time geometry itself.HI In arbitrary space-times, the high 
frequency limit is completely governed by the eikonal equation, 



where x a = (x\t) are local space-time coordinates, g ab (x c ) is a given space- 
time metric and S = S(x a ). The level surfaces of S, i.e., S(x a ) = constant, 
define the characteristic or null three-surfaces (or what Arnold calls the "big 
wave fronts"!) and the S(x\t) = constant for constant t define the two di- 
mensional "small" wavefront in the three dimensional space of t = constant. 
The vector p a = g ab dbS is tangent to the null geodesies that rule the char- 
acteristic surface. Though we are basically interested in Eq.([l|) for arbitrary 
space-times, here we will confine ourselves to a study of its solutions only in 
flat (Minkowski) space. (A future paper, in preparation, will generalize the 
present material to curved space-times.) Eq.(|l|) then becomes 



The level surfaces of the solutions to Eqs.(|T|) or (||), the null surfaces, can 
be viewed as being generated by the evolution of 2-dimensional wave fronts. 
Specifically, a wavefront evolves by following light-rays that are normal to 
it, generating the null three-surface. A smooth wavefront in three-space, in 



g ab d a S b S = 



(1) 



n ab d a Sd b S = {d t Sf - {d x Sf - {d y Sf - {d z Sf = 0. 



(2) 
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general progresses, into a singular one, either to the past or the future, i.e., a 
generic null surface in space-time has singularities. The singular wavefronts 
are 2-surfaces that are continuous with existing first derivative, but where 
(piece-wise) the second derivative are singular, being either undefined or 
infinite. The structure of the singularities are generically cusp ridges and 
swallowtails. There are unstable exceptions. 

A textbook example!! of flat space singular wavefronts (and associated 
big wavefront) are from imploding triaxial ellipsoids, where an initially ellip- 
soidal wavefront is evolved inwardly, it self-intersects for some finite period of 
time and eventually expands out to infinity, becoming spherical in the limit. 

The singularities of wavefronts are also interpreted as the location of fo- 
cusing regions, where the intensity of light becomes very high. At the focusing 
regions, neighboring null geodesies meet, and the cross sectional area of the 
bundle of light-rays collapses to zero, which leads to the increase in intensity. 
Spherical wavefronts focus at a single point, (which are unstable under small 
perturbations of the front) whereas generic wavefronts trace spatial curves of 
focusing points (cusp ridges and swallowtails). 

In Section II we will review a beautiful method for giving the general 
solution to the flat-space eikonal; it allows for solutions either from arbitrary 
Cauchy data or for stationary solutions that arise from the ansatz, S = 
t - S (x,y,z). 

In Sec. Ill, we will apply the method of Sec. II to obtain global asymp- 
totically spherical null surfaces that are associated with shearing ("bad") cuts 
of null infinity, 0.0 They will be defined from the normal rays to a "bad" cut. 
This construction can be thought of as beginning with a deformed, initial, 
two-sphere in a finite region of space-time, then construct the future outward 
directed null normals to the two surface which generates a null surface and 
finally "slide" the initial two-surface along the null geodesies that generate 
the null surface, to future null infinity. This limit is the "bad" cut of null 
infinity. 

In Sec. IV, we will study the caustics and singularities of these charac- 
teristic surfaces and their associated wavefronts. 

In Sec. V, we treat the same problems of the singularities of these surfaces 
but now from an alternate point of view, namely from Arnold's theory of 
generating families.^ This treatment allows one to handle (parametrically) 
the regions of self-intersection and non-smoothness of the null surfaces. 

In Sec. VI, we discuss a generalization of the ideas presented to this 
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point. Though this generalization is primarily intended for use in non-flat 
Lorentzian space-times, nevertheless we believe that it is quite useful to see it 
in the simpler case of flat space-time; it allows for the clarification of certain 
points that would be difficult in more general situations. Specifically we will 
consider solutions of the eikonal equation that depend on two parameters - 
that are different from the two parameter family of plane waves. We will see, 
the slightly surprising result, that the singularity structure of the individual 
characteristic surfaces can be studied via the parameter behavior of nearby 
solutions. More precisely, if the two parameter set of solutions is given by 
Z(x a , /j,, ~p), the singularities of the level surfaces of Z for fixed values of (//, ~p) 
can be studied and expressed in terms of (/i, fx) derivatives of Z. These results 
become important in asymptotically flat space-times where the Z(x a ,fi,]l) 
can be chosen to represent the family of past null cones from all the points 
of future null infinity. 



II. Solutions of the Eikonal Equation 

We review a powerful method for solving the flat-space eikonal equation 
with arbitrary given Cauchy data. We begin with a solution S* of the eikonal 
equation that depends on three arbitrary parameters, i.e., 



S* = S* (x\ t, on) = x^i - ty/E(ati) 2 



(3) 



called a complete integral. A "general integral" (which involves an arbitrary 
function) can be constructed from the complete integral in the following 
manner: we first add to it an arbitrary function of the three ctj , i.e., we 
consider 



S** = S*(x i ,t,a i )-F(a i ), (4) 
with the weak condition that (aside from lower dimensional regions) 

d 2 S** 



daidcuj 



^0. 



(5) 



We next demand that dS**/dai = dS*/dai — dF/dai = 0, which implies 
that there are three functions of the form ai = Ai(x\t). (In general these 
solutions are not unique and they must be expressed on different sheets. See 
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Sec. IV for a complete discussion of this issue.) Finally, via = Ai(x\t), 
the ctj are eliminated in the S** yielding (perhaps multivalued) 

Six*, t) = S*(x\ t, A t (x\ t)) - F(A t (x\ t)). (6) 

The level surfaces of this S might self-intersect and be only piece-wise 
differentiable. 

It is not difficult to show that the S, so constructed satisfies the eikonal 
equationl! This follows immediately from the fact that 
d a S = d a S* +(d Ai S*- d Ai F)d a Ai = d a S*. 

This solution now depends on an arbitrary function of three variables, 
namely the F. The task is now to determine the F(a,i) in terms of (appro- 
priate) Cauchy data, Sc( x l ). This is accomplished as follows; define an = 
d Sc/d x l and algebraically invert it in the form of the three equations x % = 
X l (a.i). At t — to we have that 

S(x\t ) = S*{x\t Q ,A i {x\t Q )) -FiMx^to)). (7) 

Replacing all the Ai by a, and all the x l by X l (a>i), we have that 

F(a i ) = S*(X i ,t ,a i )-S c (X i ), (8) 

i.e., the free F(oii) is now expressed in terms of the free Cauchy data, Sc(x l ). 

This allows us to find (in principle - modulo algebraic inversions or 
parametrizations) solutions of the flat-space eikonal equation with arbitrary 
Cauchy data. 

There exists a special class of solutions that are not studied or found 
easily via the Cauchy problem, namely the "stationary" solutions which have 
the form 

S = t-S (x i ). 

To generate solutions of this form we modify the complete integral, Eq.|3|, 
making it a function of only two free parameters by imposing the condition 
that S(«j) 2 = constant; (for convenience chosen as S(«j) 2 = 1/V2). 
We then write the modified complete integral as 

S*=x a £ a (CX) (9) 

where 
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4(C, = ^ [(i + CO, -(C + 0, -»(C - 0, (i - CO] (io) 
V4 1 + CO 

The complex number £, which plays the role of two of the independent 
parameters among the three can be thought of as the complex stereo- 
graphic coordinate on the sphere; the £ a (CX) is a Minkowski null vector, 
i] ab ^ah = 0, that spans the entire lightcone as ( ranges over the sphere. EqJ| 
represents a spheres worth of different families of plane waves parametrized 
by the direction £. 

If we now take 

S** = x a £ a ((,()+a((,() (11) 
and construct dS** = ~&S** = 0, i.e., 



co = x a m a (C,C) + 5«(C,C) = 0, (12) 
ZJ ee x a m a ((,() + M((X) = (13) 

where 



m a (C,C) = 5£ a (C,C) = (l + CC)^ 1 ^ (14) 

- - - - BP (C 7) 

m a (C,C) = 5£ a (C,C) = (l + CC)^ 1 ^. (15) 

For generic a(C,0? Eq-QTID can be solved for 

C = T(x,y,z), (16) 

where again these solutions need not be unique and must be expressed on 
different sheets. (See Sec. IV for a full treatment of this problem.) Note that 
Eqs. (|i~2|) do not depend on the time t and hence T is a function only of the 
spatial coordinates. When the T(x,y,z) is substituted into Eq. (|Tl|) , i.e., 

S(t, x ,y,z) = -j=- S (x, y> z ) = -L + x%(T, T) + a(T, T), (17) 
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we have a solution of the eikonal equation depending on an arbitrary function 
of two variables, a(T, T). The level surfaces of S could in general self-intersect 
and be only piecewise differentiable. 

The procedure of beginning with a complete solution and obtaining the 
general solution via the two (or three) variable arbitrary function is geomet- 
rically equivalent to the construction of an envelope from the family of plane 
waves as the two (or three) constants in the complete solutions are varied.lll 
Since in this work, we will only be interested in individual null sur- 
faces and their properties, we can and will confine ourselves to the level 
surfaces of the solutions of the form given in Eq. (|T7|) . 

III. Null Surfaces Generated by Normals to Two-Surfaces 



We want to give a slightly different geometric interpretation to the method 
of the previous section for generating the stationary solutions of the eikonal 
equation. Given any (spatial) two-surface (for example consider any two di- 
mensional slice of the past light cone of an arbitrary space-time point), the 
normal rays to the surface (either the outgoing or incoming ones) generate 
a null surface. In this section we will consider a particular case of this con- 
struction where this past lightcone is taken to be the future null infinity, 3 + , 
of Penrose.0 

The future null boundary, 3 + , (the endpoints of future directed null 
geodesies) of any asymptotically simple space-time is a null surface with 
topology R x S 2 . A choice of Bondi coordinates (u, £, C) can De assigned to 
3 + , where ( = cot(9/2)e l< ^ for the S 2 sector. The intersection of the fu- 
ture lightcone <£ x of a point x a with 3 + is a two-surface, locally imbedded in 
R x S 2 ; it can generically be described locally by u = Z(x a , (,(). The two- 
surface is referred to as a lightcone cut,0 whereas the function Z(x a ,(X) 1S 
referred to as a lightcone cut function, and is a two-point real function on 
the space-time and the boundary, 3 + . 

Remark 1 Though for this work it is irrelevant, we mention that the light- 
cone cut function Z(x a , £, Q is one of two fundamental variables in a refor- 
mulation of general relativity via null surfaces^ It encodes all the conformal 
information of the space-time. . 
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In the remainder of this work we will confine ourselves to flat space-time 
where (modulo Poincare transformations) a natural choice of Bondi coordi- 
nates (u n , C, C) exists; the u n = constant is constructed from the intersection 
of the future lightcone, <£(t, 0)0)0 ), of the spatial origin, at time t = u n , with 
3 + ; the (C, C) are j us t the null directions, at the origin, carried along by 
the null generators of the lightcone. Using Cartesian coordinates x a for the 
space-time and these natural Bondi coordinates,!^ the lightcone cuts can be 
described as 

u n = x a £ a ((,() (18) 

where £ a represents the covariant version of a null vector £ a with Cartesian 
components given as 

nc,0 = * ((i + CC),(C + C)^(C-C),(CC-i)) (19) 
V2(i + CO 

= — — (1, sin^ cos0, sin 6* sin 0, cos 9) . 

Adding the radial coordinate r, this natural choice of Bondi coordinates 
is identical to the standard null polar coordinates (u n ,r, (,() given by 



x 



ut a + r£ a ((,(), t a = 72(1,0,0,0) (20) 



Note that Eqs.fllJ) and (fj) are identical with Eqs.(| ) and (0) though 
their meanings are different; Eq.(|l8|) has the dual meaning of being, for fixed 
value of the x a , the lightcone cut of 3 + and also, for fixed values of (u n , (, (), 
it describes the plane wave (null surface) intersecting the time axis at t = u n 
in the direction of (CO- 

By Eq. (|l8l) , the lightcone cuts of any points x a = (t, 0,0,0) along the 
time axis take, as we mentioned earlier, a constant value on 3 + , namely, they 
are the constant-w n slices. The natural Bondi cuts are lightcone cuts as well. 
By following inwardly the null geodesies that leave the natural Bondi cuts 
orthogonally, we find no focusing other than at the apex (on the time axis) 
of the lightcone. 

By a slight modification of the above we can find other null surfaces 
leaving 3 + that have much more complicated focusing properties than that 
of a lightcone. If we consider the one parameter family of cuts of 3 + given, 
say by, 

u n = -a((,()+u (21) 
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where a(C, C) is a given but arbitrary regular function on S 2 , and u is a 
parameter on R, we can ask for the null surfaces generated by the null normals 
to the family of cuts. 



• Note that Eq. (|2l|) can be rewritten as u = u n +a((, () and reinterpreted 
as a (Bondi-Metzner-Sachs) supertranslation0 between the coordinates 
u n and ii on J + . 

We now construct the null surface formed by the normal rays to the cuts, 
Eq.(|2~l]), determined by u = const.; replacing the u n in Eq.(pT|) by the null 
planes Eq.(|l8|), \u n = x a £ a ((,Q], we have 

u = x a e a (U) + a((,() (22) 



which is identical to Eq. ([TTD . The envelope formed from all the null planes 
that are normal to the family of cuts are found by setting to zero the ( and 
C derivatives of Eq.([22|) and eliminating the ((,() from (|22"D, a procedure 



identical to that followed in the previous section to obtain Eq . (|17|) , i.e., we 
now have the one parameter family of null surfaces 

S** = u = —j= — x7i(T, T) + a(T, T) (23) 
v2 



with C = Y(x,y,z), a solution of Eqs. (|12]) , (JT3|) . The procedure of setting to 
zero the u and uJ in 



u = x a m a ((,() + <ky((,() = 0, (24) 
ZJ = x a m a (C,C) + 5«(C,C) = (25) 

selects the null ray at each point of 3 + that is orthogonal to the cut given by 
E q .(0). 

• Note that a can be chosen to contain only spherical harmonics of order 
I ^ 2 since any I = 0, 1 components of a could be absorbed by x a £ a 
with no modification other than displacing the origin of the coordinates 
x a , since £ a is precisely the collection of spherical harmonics of order 
and 1. 
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We can give a parametric description of the family of null surfaces, 
Eq-(H), by the following procedure: we consider the four functions 

C = T(x, y, z) x a m a ((, () + <Sa((, = (26) 
C = T{x, y, z) & x a m a (C, () + 5a(C, C) = (27) 
u = u(t, x, y, z) = {x a £ a ((, C) + a(C, OIItTt (28) 

r = r(x,y,z) = {x a (n a - () + 55a (C, C)}ItT ( 29 ) 

and consider them as a coordinate transformation between the {u, £, C) r } 
and the {x a }. We have used m a = 5£ a , rh a = M a , and n a = 55£ a + i a . From 
the fact that (£ a , m a , m a , n a ) form a null tetrad for every fixed value of (£, (), 
this coordinate transformation can be readily inverted into the form 

x a = (u — a) (n a + 1 ) + (r - 55a)£ a + (5a)m a + (5a)m a (30) 

This relationship can alternatively be looked on as the parametric version of 



the one-parameter family of null surfaces, Eq.(|23|), (where, for fixed u, the 
(r,£, C) parametrize the surface) or as the coordinate transformation between 
the {x a } and the null-geodesic coordinates, (u,(X, r )'i u labels the null sur- 
faces, the pair ((, Q labels null geodesies (via their intersection with 3 + ) and 
r is an affine parameter along the null geodesies. That this is true can be 



easily seen from the parametric form, Eq.(|30|), by simply constructing 

dx a 



f// , r(C,C) (31) 

and observing that £ a is a null tangent vector with affine normalization. 

The transformation between the x a and the (u, r, (, () breaks down when 
the Jacobian of the transformation, Eq.(|30|), vanishes, i.e., when 



m^7xx)^ r _8a8a = - (32) 

Geometrically, this is where the null surface develops singularities. In the 
projection to the three-space (x, y, z) it is a two surface; the "caustic surface". 
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To see this explicitly, we return to Eq.(|30D where we have (for fixed u) that 
(x, y, z) = x l = X l {r, C, C), i.e., are known functions of (r, CjC)- If the r in 
X 1 is replaced by the r from Eg. ([32]) we have the parametric form of the 
caustic, 

x* = X\r((,()XX)=X\(,C). (33) 

We will return to this in the next section. 

It is interesting to note that the coordinates (u, (, (, r) represent a type of 
null coordinate system that we could call asymptotic null-polar coordinates 
which are the flat space case of an interior Bondi coordinate system, @0 
i.e., the extension into the interior of the space-time of the Bondi coordinates 
(it,e, C) 011 3 + - They differ from the standard null polar coordinates by the 
fact that the null geodesies that rule these surfaces possess non-vanishing 
shear while for the standard ones the shear vanishes. 

The complex shear is defined as a = M a M b V a Lb, where Lb is tangent 
to the null geodesies and M a is complex null, orthogonal to Lb and such 
that M a M a = —1. In our case, because of Eqs.(|30|) and ( |3ll ) and the fact 
that (£ a , m a , rh a , n a ) forms a null tetrad, we have that Lb = h and M a = m a . 
Furthermore, the gradient of 4 is V a £b = m a (,b +"TaCfc an d thus a = —m b (,b ■ 
To obtain the derivative of ( along m b we take the gradients of Eqs.(|24]) and 
pal ) which yields 



m b {(, + - la) + 35«}C, fe +5 2 «C, 



r2 t 



0. 
0. 



m b(C, + {x a (n a - l a ) + SSajCfe +5 a( }b 
Using Eq.(f29|) and contracting Eqs.(|34]) and (|35|) by m b , we obtain 



(34) 
(35) 



rm b (, b +5 2 am 6 C, 6 



0. 



-1 +rm b (,b+d 2 am b (, b = 0. 



(36) 
(37) 



By eliminating the m b (,b from these equations, we find 



(r 2 -5 2 «5 a)m b (, b +tfa = 



or 



a 



a 



r 2 - a°a° ' 



(38) 
(39) 
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where cr c 



5 2 ct. This is also a confirmation of the Sachs theorem on the 



transformation of the asymptotic shear, a , under a BMS transformation. 
Eq.(|3~9|) represents a special (non-twistinglcase of a more general result valid 
for generic null congruences in flat spaceHI 



We point out that the flat-space line element, using Eq. 
be expressed in terms of these shearing Bondi coordinates a 



can easily 



ds 2 = r] ab dx a dx b 



2 _ d C _ d( 



2du ^du + dr — 55 a— — 55 a ^ 
- ^(5 2 a d( 2 + ta d( 2 ) - 2(r 2 + 5 2 «5a)^, (40) 



where P = 1 + CC ■ The determinant of g a b is given by \g\ = -pz{r 2 — 



5 2 a 5 a} 2 = whose vanishing agrees with the vanishing of the Jacobian, 



Note that the asymptotic r = const. => oo sections, at u = const., become 
metric spheres. 



• It is perhaps interesting to speculate on the use of Eq. fllPP as the 
Minkowski space lowest order term, in perturbation calculations, for 
solutions of the Einstein equations. 

We complete this section by showing how a null surface can be constructed 
explicitly from the normals to an arbitrary space-like two surface, &, in a 
manner virtually identical to those constructed from a cut or slice of null 
infinity. 

We begin from Eq.(^) 

u = o = z%(c,C) + «(C,C) 



with Eqs. (|2J and [25]). The issue is, given the surface &, how is one to 
choose a((, £)? 

First (3 is defined parametrically, x a = Xq(C, C), where the parameters are 
chosen as follows: consider a time-like world line at the spacial origin and 
the family of light-cones centered on the line. The null geodesies ruling these 
cones are labeled by their directions (£, () on the sphere and coincide with 
labeling of the generators of null infinity. The points on & are now (locally) 
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parametrized by the labels of the null geodesies passing thru those points. 
With this parametrization the function a((, () is defined by 

a(c,C) = -*8(C,C)4(C,C) 

so that Eq. ([23|) becomes 

u = = (x a -x a ((,())£ a (U) 
and |25|) become 

(x a - x a (C, 0)m a ((, C) - 4(C, OK = and c.c. (41) 

We see that the null surface so defined goes thru (5, i.e., thru x a = Xq((, (). 
To see that it is also normal to 6, we notice that at & the first two terms of 
Eq.(|4]) cancel out and we are left at & with 

Thus as was claimed, the tangent vectors to (5, namely 5xg(C,C), are 
normal to the null tangent vectors to the null surface, £ a . 

We see that the earlier construction of null surfaces from cuts of null 
infinity actually includes those constructed from finite surfaces. 

IV. Wavefront evolution and singularities 

In the previous section, we mentioned that the null coordinate system 
broke down and the associated shearing null surfaces developed caustics at 
the points where r 2 = 5 2 a 5 a. 

Here we focus our attention on the two dimensional wavefronts associated 
with these null surfaces. We show that the wavefronts develop singularities, 
and we locate the singularities via the standard method of singularity theory, 
and via our lightcone cut approach. The evolution of these singularities as 
the wave fronts evolve become the caustics. 

A wavefront is, by definition, the intersection of our null surface, u = w , 
with a constant-time t surface. This represents an instant in the progression 
of a wave. In our case, this requires fixing the time coordinate x° = t in 
Eq.(|30|), and solving for 

r = V2t - 2u + 2a + 55a. (42) 



and Eqs.(24 
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The remaining coordinates (x, y, z), using Eq.(|42|) to eliminate r, trace a two- 
surface (a "small" wave front) in the Euclidean three-space, parametrized by 
(CC) [ or (0)0) under the transformation ( = e 1 ^ cot 0/2]: 



1 

x = —= 
V2 

i 

v = 71 

1 

z = —r= 

The map (C 



(V2to - 2u + 2a) 5 + £ + Sc/ J ^ 



;i + CO 



i + CO 

2x 



+ 5a- 



(v^o - 2m + 2a) 5 - 5a ^ " 



2t - 2m + 2a) 



(1 + CO 

(CC - 1) 



+ 5a- 



i + CC) 

2C 



+ 5a 



5a- 



;i+co 

:i+c 2 ) 
;i + co 

2C ' 



(43) 
(44) 
(45) 



1 + CC) (1 + CC) (l + COJ 
x, y, z) is singular at points where the Jacobian matrix 



5x 5x 
by by 

bz 5,2 



(46) 



drops rank, from 2 to 1 or 0. The drop in rank takes place if the three 
2-determinants vanish simultaneously: 

(47) 
(48) 
(49) 



5x5y — 


5y5x 


= o, 


bybz — 


bzhy 


= o, 


5z5x — 


5x5z 


= 0. 



Since 



bx a \ Uit 
bx a L t 



b 2 am a + rm a , 

<2 



5 am a + rrh a , 

the explicit expressions of the 2-determinants are as follows 

5x5y — bybx = 



r 2 - b 2 ab 2 a){m x m v - m x m y ) 



bybz — bzby = (r 2 — 5 2 a5 a){m y m z — m v m z ), 
bzbx — bxbz = (r 2 — 5 2 a5 2 a)(m 2 m x — rh z m x ). 

Thus, all three determinants vanish at points where 



D 



5 2 a5 2 a 







(50) 
(51) 

(52) 
(53) 
(54) 

(55) 
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which, with r from Eq.fl42|), determines a curve, the wavefront singularities. 
Note that the evolution of the wavefront singularities (obtained by varying 
t ) yields the caustic surface. 

Since a is a regular function on the sphere, so is 5 2 a 5 a; therefore, 
Eq.(|55|) (with r given by Eq.([42|)) admits solutions ( only within a finite 
interval of time t. Thus the wavefronts are singular only during a closed 
interval of time. On the other hand, at very long times the wavefronts become 
spherical, which follows from the line element Eq.(^0|). 

The wavefront singularities (curves) are places where neighboring null 
geodesies meet. We have a null surface x a (r, (,() foliated by null geodesies. 
At every fixed value of r, there are two connecting vectors 5x a | u r and 5x a | u r . 
The null geodesies in this congruence meet wherever the area orthogonal to 
the congruence, spanned by the connecting vectors, vanishes. The connecting 
vectors are, explicitly, 

d x a \ ur = - 55 2 a r + a 2 am fl + rm a , (56) 

5a; a | Uir = - dtat + rm a +tam a , (57) 

The area spanned by the connecting vectors (calculated from their skew 
product) is simply D = r 2 — 5 2 a5 a. The vanishing of this area takes place 
at exactly the points given by Eq. (|55|). 
We close this section with two examples. 

Example 1: a = Y20 = 3 cos 2 9 — 1. Due to axial symmetry, the wave- 
fronts and their singularities for this choice of a can be completely worked 
out analytically, which gives insights into more general cases. The wavefronts 



at a given time t are given by 

x = —j= sin 6 cos 0(v / 2t-2M-2-6cos 2 ^), (58) 

y = —= sin 9 sin 0(v / 2t-2M-2-6cos 2 ^), (59) 
V2 

z = -^cos9(V2t-2u+10-6cos 2 6). (60) 



These are axially symmetric. For a closed interval of time, all the wave- 
fronts are singular. For early and late times, however, the wavefronts are 
smooth. 
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The singular points are located by eliminating r from Eqs.([55|) and fl42"D , 
yielding in this case, the two solutions or "sheets" 



There is a solution 9 only at times y2(u — 5) < t < \/2(4 + u). This 
is the interval where every wavefront is singular. A smooth wavefront and 
its corresponding profile at a late time are shown in Fig. 1. A wavefront at 
a time when both the cusp ridge singularities and the z -axis singularities 
are occurring, and its corresponding profile, are shown in Fig. 2. In Fig. 3, 
we have a later wavefront and profile with only the cusp ridge singularity. 
In Fig. 4, we display the evolution of the singularities forming the caustic 
2-surface. 

The high symmetry of this case is responsible for the lack of resemblance 
of the singular points on the z-axis with standard cusps. At these points, null 
geodesies labeled by different, but neighboring values of 0, meet. This is clear 
from the fact that dx/dcj) = dx/dcj) = dz/dep = at these points, therefore 
the vector that connects geodesies with different values of ep vanishes. 

In order to make a comparison, Fig. 5 shows a wavefront in the evolution 
of an imploding ellipsoid of revolution which is very similar to that of Exam- 
ple 1. In this case, an ellipsoid of revolution sends an incoming wavefront, 
which develops singularities during a certain interval of time. The standard 
cuspoidal ridges are clearly visible as rings at both ends of the figure. How- 
ever, the crossover points in between are also singular, of the same type of 
singularity as that one developed in our example. Assuming a speed of light 
of 1, the formulas for the imploding wavefront in this case are 



(y/2t - 2m + 10 - 18 cos 2 9){y/2t - 2u - 2 - 6 cos 2 9) = 0. 



(61) 



x 



y 



z 




(62) 



(63) 



(64) 



Example 2: a = Real(Y 2 i) = (C + CXCC - 1)/(1 + CC) 2 - In this 
case, there is no advantage in writing the wavefronts explicitly. 
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However, they can be plotted with ease, displaying the typical 
singularities of three-dimensional wavefronts, namely swallowtails 
and cusp ridges. Cusp ridges are clearly visible in Fig. 6, which 
represents a wavefront at u = and t = 1.5. Swallow tails are 
exemplified in Fig. 7, which represents another wavefront in the 
same u = family, but at a later time of t = 2.35. Locally all 
swallowtail have the form of Fig. 8. 

Both Figs. 6 and 7 compare remarkably well with a typical 
wavefront in the evolution of a triaxial ellipsoid, in which case 
an ellipsoid (x/a) 2 + (y/b) 2 + (z/c) 2 = 1 emits a wavefront of 
light inwardly, which develops singularities for a period of time. 
A typical singular imploding wavefront is shown in Fig. 9. The 
formulas for the imploding triaxial-ellipsoidal wavefront are the 
following: 



x = a sin 6 cos 0(1 — 
y = b sin 9 sin 0(1 — 



t 



! a/ (sin 8 cos (J)/ a 2 ) 2 + (sin 6 sin 4>/b 2 ) 2 + (cos8/c 2 ) 2 ' 

t 



6 2 A/(sin 6 cos <f)/a 2 ) 2 + (sin 6 sin <p/b 2 ) 2 + (cos 6*/c 2 ) 2 ' 

z = ccosfl(l- = ) (65) 

Y (sin 6 cos0/a 2 ) 2 + (sin# sin <p/b 2 ) 2 + (cos8/c 2 ^° 



V. Generating Families 

In this section we will study the subject of the caustics of the null surfaces 
and the wavefront singularities via an alternative method, namely from the 
use of generating families for the construction of Lagrangian and Legendre 
submanifolds (developed by V. I. Arnold and his colleaguesfilSB'Q) associated 
with cotangent and contact bundles over space-time. The value of this treat- 
ment is that it allows one to deal (via a parametric representation) with the 
regions of self-intersection and non-differentiability of the null surfaces. 

We first give a brief review of a special case of this theory that is adapted 
to the problem of null surfaces in four-dimensional space-time. Consider a 
four dimensional Lorentzian manifold, (with local coordinates, {t, x 1 }) foli- 
ated by the constant t surfaces. Now consider the x % as the coordinates of a 
configuration space M and pi as the conjugate momentum so that we have 
the six dimensional cotangent bundle T*M, with local coordinates (x\pi). 
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We now describe the construction of a three dimensional submanifold of T*M 
(a Lagrangian submanifold - a maximal submanifold such that the symplectic 
form, restricted to it, vanishes) that plays a fundamental role in the discus- 
sion of the singularities of wavefronts and their associated caustics. We begin 
with a general description without any particular choice of dynamics, later 
restricting ourselves to null geodesic motion. 

Choose a scalar function (determined later from the dynamics), re- 
ferred to as a generating family, of the six variables 

f = F(t,x\ c,C); (66) 

the x % are spatial points, the ((, Q are parametric labels (for convenience we 
are using a complex representation) for points on a given spatial two-surface, 
5, i.e., we have a two-point function, while t is the time for the (dynamic) 
particle to go from a point on s to the point x\ For a constant value of /, 
we consider Eq.flBlf) as defining t implicitly as a function of (x\ £, Q, i.e., 

t = T(f;x\U) 

or simply 



t = T(x\CO- (67) 

Note that T might be a multivalued function of its arguments, in which case 
it must be considered separately on the different sheets. 

We now ask for the relationship between the (x\(,() when T is an ex- 
tremal under variations of the (C>C)> i- e - 5 we require that 

OT/dC = OT/dC = (68) 

which in turn forces 

dF/dC = dF/dl = 0. (69) 

Finally a rank condition is imposed on the choice of F; the following 2x5 
matrix must have rank two 



" d 2 F 


d 2 F 


d 2 F I 


d*F 


acac 


<9f dx' 


d 2 F 


d 2 F 


d(dC 




d^dx' 



(70a) 



The meaning of this condition is that the two equations (|69] or |68|) can be 
solved (locally) in at least one of a variety of possible ways for two of the five 
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variables {x\ C,C)j often it is necessary to solve them in the different ways in 
different regions. We then have three different possible cases: 

1. ( = T(a;*), ( = T(x l ); the simplest of the three cases. It allows F to be 
treated as a function of just x l . In the other cases F must be treated 
parametrically. 

2. C — iC)i % J — X J (x A ,(), where x A are any two of the three x % 
and x J is the third one - or the conjugate version, ( = ^l(x A , (), x J = 
X J (x A X). 

3. x A = X A (x J , C, (), where again x A are any two of the three x l and x J 
is the third one. 



Case I can occur when the determinant 



Q2p Q2p^ 

D= ~% P§ t^O; (71) 

D — when there is a critical point of the Lagrangian map defined shortly. 
The Lagrangian submanifold obtained from F is defined in the following 

way: 

First we have Pi = (Note that this involves only the explicit x % de- 
pendence in F since any implicit dependence, via the (C,C)> does not enter 
into the definition of Pi because of Eq.([6^).) Now depending on which case, 
#1, 2 or 3 is relevant, we eliminate two of the five variables (x\ (, (), in the 
Pi = j^j(x\(,(). This leaves the result that the six coordinates of T*M, can 
be expressed in terms of three parameters, thus defining a three dimensional 
submanifold in each of the cases; 

• In case 1, Pi = Pi(x l ), x l = x % \ the three parameters {x 1 } = \ a - 

• In case 2, pi = Pi(x A ,(), x A = x A , x J = X J (x A ,Q; the three 
parameters {x A ,(} = x a - 

• In case 3, pi = Pi(x J , CjC)j x — X A (x J ',(,(), x J = x J ; the three 
parameters {x J ,(X} — X a - 
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To simplify the discussion, we have referred to the three parameters, in 
each of the cases, simply as x a - I n eac h case we thus have 



x l = X\ X a ), Pi = Pi(x a )- (72) 

Of immediate relevance to us is the projection (Lagrange map) of the 
Lagrange submanifold into the configuration space which becomes in each of 
the cases, x l = X l (x a ) or 

• In case 1, x l = x\ trivial diffeomorphism 

• In case 2, x A = x A , x J = X J (x A , () 

• In case 3, x A = X A (x J , (, £), x J = x J 



The caustics of this problem are the regions in the configuration space 
where the mappings #1, 2 or 3 have rank two or one; i.e., where the Jacobian 
of the mapping vanishes. They occur when the determinant D — 0. The 
inverse image to the caustics in the parameter space are referred to as the 
critical points of the Lagrange map. It is clear that in case 1, the Jacobian 
is one and rank reduction can only occur in cases 2 and 3. 

This treatment of caustics can be extended into the full four-space by 
eliminating, in the expression for t = T(f; x' 1 , £, £), or m the implicit version, 
/ = F(t, x\ C, C), two of the five variables (x\ £, () via the cases #1, 2, 3. This 
results in t now being a function of the three parameters, x a - Though we 
will not need the full theory here, this construction leads to a seven dimen- 
sional manifold, (t,x\pi) {an example of a contact manifold} and a three- 
dimensional submanifold of the contact manifold (a Legendre submanifold) 
defined by 

t = T( X a ), x* = X\x a )i Pi = Pi(x a ) (73) 

as well as the Legendre mapping from the Legendre submanifold, to space- 
time, (t, x % ) 

{t,x\ Pl }(x Q )^t = T( X a ), x l = X\ X a ) (74) 

a three-surface in space-time - the "big wavefronts" in Arnold 's language - in 
our case a null surface. The singularities of the map Eq. ([74|) , where the rank 
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drops below three, are the null surface singularities. These singularities, at 
fixed t, are the ("small") wavefront singularities of the previous section. 

We now return to the question of the determination of the function / = 
F(t, x\ C, C) °f Eq.fl66|) for use in the study of null surfaces. Our choice will 
be, from Eq.(pTT|), 

F(t,x\(,() = S**(t,x\U) ee x a l a ((X) + a((X). (75) 
There are three independent reasons for this choice; 

1. It was the method of generating an arbitrary null surface from the com- 

plete solution, x a £ a ((X)i see Sec. II.. 

2. It was the method for the construction of a null surface such that the 

generators were orthogonal to a given two-surface; see Sec. III.. 

3. It arises from a variant of Fermat's Principle of stationary time: Consider 

a timelike worldline, £, (in a Lorentzian space-time) of the form, in local 
coordinates, (x l = constant, t varies) and a two-dimensional space- 
like surface, s((,Q- Assume, locally, that from every point of s((,() 
there is a null geodesic that reaches £ at a time t = T(x\(,(), then 
t is extremized by those curves that are normal to s(C,C)- This result 
(which will be described in detail elsewhere) follows from Schrodinger's 
derivationll! of the gravitational frequency shift. 

Note that the rank condition on the matrix, Eq. ( |7Ua| ) is satisfied by direct 
calculation. 

From the discussion of generating families, we see that the treatment 
of the null surfaces that we proposed, in Sees. II and III, namely to solve for 
the C = T(x, y,z), was really only valid for Case 1, but we actually used a 
version of Case 3 where the additional parameter r was introduced in order 
not to single out any particular cartesian coordinate. Case 1 broke down 
precisely on the caustic, given by 

r 2 = tfata (76) 

which is where Cases 2 and 3 must be applied. 

What follows is a straightforward application of Case 3 of these ideas - 
on one patch. For completeness we repeat some of the earlier steps 



22 



Starting with 

u = S**(x a , C, = x a UC + «(C, C) (77) 



then 



5S** = -=^- — ^[CV- W-2zQ +5a, 

55** = - 7 =-L_[C 2 iy-W-2zC]+3a, (78) 

with PF = (x + iy). From 3S"** = and 5S 1 ** = 0, we obtain, (from Case 3, 
where x A are (x,y) or (W, W)) that 

w = ^;s^-^s > (79) 
w = 



From Eq.([77|) we have that 

= 

Py = 





+ C 2 3«) - 2z( 




1-CC 




+ C 2 9a) - 2z( 




i-CC 


dS** 


c + c 


dx 


v^(i + CO ' 


dS** 


i(C-C) 


dy 


v^a+co' 


dS** 


(i - CO 



P* = = r- - (80) 

Taking p = V2(p x + , we obtain 

P= ^77- (81) 

l + CC 

The Lagrange submanifold, parametrized by (z, CO; i s gi ye by 

2 = 2, (82) 

= v^(5a + c 2 a a) -2,c (g3) 
i-CC ' 
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W = V2(8a + r9«)-2< (g4) 

p - -T^C- (») 

P = (86) 

ft = < 1 - Cfl - . (87) 



The projection to the configuration space is given by 

z — z, 

V2(M + ( 2 M) - 2z( 



W = 



i-CC 



- V2(da + C 2 5a) - 2z( 

W = i s ; -. (88) 

The critical points of the Lagrange map are obtained from the condition 
that the Jacobian 

d(z,C,0 

vanishes. The vanishing of J is equivalent to D = r 2 — 5 2 a 5 2 a = 0. 

To construct the Legendrian submanifold (in the seven dimensional con- 
tact space, (t,x\pi)) we take the generating family u = S**(x a ,(X) where 
u is constant and solve for t expressing the contact coordinate t in terms of 
the three parameters (z, (, () by 

t = (1 + CC) + V2[u-a{C,Q\ (90) 

with 



W 
W 



V2(ba + ( 2 M) - 2z( 
V2(Sa + C 2 9a) - 2z( 



(91) 
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The full Legendre submanifold is then given by Eqs.(|2|), (g3|), flS§, (||), 
©, @ and ©. 

Note that this entire construction, using Case #3, was valid where 1 — 
CC 7^ (or equivalently where p z ^ 0). To include the region where p z = 0, a 
different choice of parametrization would be necessary, e.g. (x, (,(), which 
is valid in the region where p x ^ or (y, (, Q, valid where p y ^ 0. 

Using the example #2, from Sec. IV, given by 

S**(x a , C, C) = x a l a (t, C) + a(C, C), a = - (1 ~ C Plit^ ( 92 ) 

(i + CC) 

with 

55** = o {x(-i + f) + iy(i + c 2 ) - 2z(] (1 + CO = V^[l + C 3 C - 3C(C + 0] 

(93) 

and 

5 = & {x(-l + C 2 ) - iy(l + C 2 ) - 2<] (1 + CO = V2[l + C 3 C - 3C(C + C)] 

(94) 

one could try to solve for different pairs from the set (x,y,z,(,0- When 
D ^ one could always solve for (CO 5 though in general there would be 
more than one solution; i.e., for fixed (x,y,z) there would in general be 
more than one ray going thru that space-point, either at the same or at 
different times. Alternately one could try to solve in different regions for 
(x,y), (y,z), (z,x), etc. Solving for (x, y) or W we have that the Lagrange 
map (from the (2, C, C) parameter space) becomes 



w = v / 2(2C 2 +4CC-l-C 2 C 2 )-2zC(l + CO m 

i-C 2 C 2 

W = v / 2( 2 C 2 +4CC-i-C 2 C 2 )-2<(i + CO (96) 

i-c 2 c 2 

which in turn becomes the Legendre map when the contact coordinate t is 
added in; 



t = ^{u+ (1 _ ^ - ^ )2 [4CC(C+0-^(l+3CC + 3C 2 C 2 + C 3 C 3 )]}- (97) 
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Though none of the details of this analysis is particularly enlightening, it 
nevertheless shows how in principle one constructs the Lagrange submanifold 
and the Lagrange map even in the presence of the caustics. 

VI. Families of Foliations 

In this section we will generalize, in the following sense, the ideas of 
Section II. Recently there has been a reformulation of General Relativity, 
referred to as the Null Surface Formulation (NSF) where the basic idea has 
been to use a family (a spheres worth) of null foliations of space-time, so that 
there are a spheres worth of null surfaces passing thru each point of space- 
time.. These surfaces are described as the level surfaces of the function 

u = z{x a XX); 

x a are the local space-time coordinates and (£, Q are the complex stereo- 
graphic coordinate on the sphere which labels the family of foliations. The 
function Z, for every fixed value of (£, £), satisfies the eikonal equation, 

g ab d a Zd b Z = 0. (98) 

Knowing these families of foliations one can construct the (conformal) 
metric in terms of Z. The idea was then to express the Einstein equations in 
terms of these surfaces, i.e., in terms of Z and a conformal factor. Though this 
was successfully accomplished, a technical difficulty in fully understanding 
the equations arose due to the fact that the null surfaces developed singu- 
larities (caustics) and self-intersections. It was clear that the development of 
caustics was a generic feature of the equations but it was not at all clear how 
to see and study their existence directly in terms of the function Z(x a , (, () 
and its derivatives. In this section, we will study, in flat space, the construc- 
tion of such families and show explicitly how to calculate the structure of the 
null surface singularities (the caustics and wave-front singularities) directly 
in terms of the Z function. 

Locally (up to first derivatives) there is no direct curvature involvement 
in the eikonal equation, so that the form of the caustics in terms of Z should 
apply equally in curved space as in Minkowski space. The results obtained 
here for Minkowski space will thus likely apply to the curved space situation. 

Starting with the two-parameter family of plane waves used earlier, Z (x a , (, 
x a £ a ((, £), we will first construct a general two-parameter family of solutions 



26 



to the flat-space eikonal equation, Z(x a , fi, fi), with (//,//) parametrizing the 
sphere; we then study the singularities and caustics of this new family. 
We begin by generalizing Eq.(p]), namely 
S** = x a £ a ((,0 + a(<;,() 
by writing a as a function on S 2 xS 2 ; i.e., as a = a(C,C HiV) an d then 
repeating the earlier procedure of setting to zero, the 5 and 5 derivatives 
with respect to the (C? C) j i- e -> considering 

u = Z**(x, C, C, m) = x a £ a (C, + «(C, C, M, V) (99) 

and 

a f z** = \z** = o (ioo) 

and then solving them (when possible) for (£, () obtaining 

C = T(x,y,z, n,~P), C = T(x,y,z,fi,Ji) (101) 

so that when substituted into Eq. (|9^) we obtain the new family which de- 
pends on the choice of a((, (, fj,, /J); 

Z(x a ,fi,]l) = x a l a (T,T) + a(T,T>,7J). (102) 

(Alternatively we could use the different cases of Sec. V, when one can 
not solve for (£, ().) It is obvious from the previous discussion that Eq. (|102|) 
satisfies the eikonal equation for each fixed value of (/z, /x). All we have done 
so far is create a new spheres worth of null foliations (wavefront families) 
of Minkowski space - different from the plane wave case of S = x a £ a ((X)- 
As in the earlier sections we could have analyzed the null surfaces for each 
value of (/i, /l) separately but now in this generalization the null surfaces are 
smoothly connected to each other through the variable (/i, fx) and it becomes 
of interest to see the development of the caustics via the variation of the 
(/i, /I), or through the (/i, Ji) derivatives. 

Remark 2 We will use, respectively, the notation (d^, d^j for the eth and 
ethbar derivatives with respect to the variables (fi,~fi) and (5^5^) for the 
variables ((, £). 
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We begin by defining several derivatives of Z\ 



to = d[j,Z, to = b^Z (103) 

R = \^Z (104) 

A = 5jZ, A = 5>. (105) 

A level surface of Z, with fixed (/i, /Z), is ruled by null geodesies, whose 
tangent vectors are given by £ a (T, Y); A pencil of rays defined from a pair 
of geodesic deviation vectors (from a given geodesic) has an area A that can 
be given@ up to a proportionality by 

A = K = (106) 



(1-AiAi) 

where K is a constant determined by the initial area and 
and 

A 1= tt~ 2 £ a A, a . 

The derivation of Eq. (|106|) is lengthy and will not be given here. It will 
however be shown, in this case to be proportional to the area. 

We now want to see the behavior of u>, A and R, as well as the area A, 
in the neighborhood of a caustic. 

By direct calculation we have, from Eqs.( |102| and |100|) , that 



uj = 5 M a, to = 3 M a, (107) 

and hence is singularity free. 

After a rather lengthy calculation, using Eqs. (|99|) and ( p.00|) , we obtain 

R = B M 3 M a + — {(9 c 9 At a)[(5 c 5 M a)(5^) - (S c o>)(9 c 9 c u)] (108) 
+ (3 ? 3 A1 a)[(3 ? 3 M a)(a 2 M) - (5 f 5 M a)(3 c 3 ( ?i)]}. 

where 
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D = (8 c 5 c u) 2 - 5^5 c u. (109) 
In a similar way we obtain 



A = a^ + -{(5 c 5 M a)[(5 c a M a)(5 cM )-(5 c a M «)(5 c 5 cM )] (110) 
+ (5 ? 5 A1 a)[(5 f 5 M a)(a^u) - (5 f 5 M a)(5 ? 5 C M)]}. 



First we see that for fixed values of (/z, /i), D from Eq.( |109|) , is the same 



as in Eq. (|32|) , namely D = r 2 — S^aS^a and hence vanishes at the caustic. 
We can now see that u is regular at the caustic while both R and A have 
singularities of the form D~ x at the caustic. 

In order to find the area A we first need R, a and A, a . After a lengthy 
calculation we obtain 



n 2 



I . v,-2 



{[(diu)(d ( u) + (5 ( 5 f w) 2 ][(5 c d M a)(5 f 5 M a) + (5 f 5 M a)(5 c 5 M a)] 



D 2 LLV^C^V^C 
,^2 



-2[(8 c 9 (U Q!)(8 c 3 iU Q!)(8 c 'u) + (5 ( 5 M a)(5 (: a M a)(5 (: M)](5 ( 5 C M)} 



fill) 



and 



2 A! = TA, C 



—{[(9^) (9^) + 



C M) 2 ](8 ( 3 M a)(3 c 3 M Q!) 



-[(8 c 8 M «) 2 (8 cM ) + (8 f 8 M «) 2 (8 2 M )](8 c 8 cM )} (112) 

Though it is not immediately obvious, from Eqs.( |lll|) and ( p.l2|) one can 
show that (1 — A 5 iA i) is proportional to D 2 or that |A i| =^> 1 + 0(D) at the 
caustic. From these results we have that 

A-KD 

(113) 



A = K 



n- 



1 - A iA, 



K 



with 



a\ 



From Eqs. (|100|) and (|101| ), we have that K = K(Y,Y) from which it 
can be shown that 



29 



i.e., K is constant along the geodesic flow. If we chose K = ±K we have 
that A = D, in agreement with the area obtained in section III from geodesic 
deviation. 

Several important observations can now be made: 

1. In our particular case of flat space, we have seen that the quantities R 
and A diverge as D~ x at the caustic. It appears virtually certain that this is 
a general result and remains true in a general curved space. 

2. As was to be expected the Area of a pencil of null geodesies vanishes at 
the caustic. This is clearly true in general and the result here is a confirmation 
that the Eq. (|106|) really is the area formula. 

3. The quantity Q (which plays a central role in the NSF version of GR) 
diverges as D^ 1 at the caustic. 

4. Though A diverges at the caustic, the absolute value of its weighted 
derivative 

|A ;1 | = \n- 2 A, a e*\ 

approaches one as 1 — 0(D). From this one sees that A, a £ a diverges as D~ 2 . 



VII. Discussion 



Our main interest in the study of wavefronts and their associated null sur- 
faces, lies in our desire to understand and describe their singularity structure 
in curved Lorentzian manifolds and in particular to find the most appropriate 
variables and representations for their analysis. Though locally the classi- 
fication of generic singularities and caustics is complete and is the same in 
both flat and curved spaces]! however, in general spaces, curvature effects are 
large and must eventually be taken into account for global questions. (For 
example, the structure of the lightcone in a curved space is very different 
from that of a lightcone in flat space.) The present work is intended to begin 
this study with the description of singular, global, asymptotically spherical, 
null surfaces in flat spaces. A follow-up second paper, will be devoted to the 
same issues as here but in asymptotically flat space-times. We will see that 
beginning with a two-parameter family of solutions of the eikonal equation - 
analogous to the plane wave solutions of flat space-time - it will be possible 
to construct any other null surface and then analyze its singularity structure. 
In particular, it is possible to construct, in terms of the two parameter fam- 
ily, the light cone of any space-time point. These insights are important for 
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applications of the null surface formulation of GR.I 
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